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ABSTRACT 



The object of this work is to formulate an analytical 
model of the heat transfer phenomena in deep submergence 
submarine shaft seals. It is expected that this will prove 
useful to design engineers in determining materials and 
arrangement specifications for these seals. 

Interfacial temperature data derived from a comprehen- 
sive experimental investigation is compared with thermal 
results computed from a lumped model of the shaft seal. 

This latter analytical data is obtained by first defining 
the model, then calculating the heat source component in- 
puts, and finally utilizing an IBM 709^ computer to aid in 
solution of the resulting mathematical entanglements. 

It was found that the most significant factor in the 
overall heat transfer scheme is the quantity of bypass cool- 
ing water flow. No bypass flow results in an extremely sharp 
temperature rise with pressure; whereas only a slight amount 
of flow will reduce the interfacial temperature rise to ac- 
ceptable levels. It was also found that with a combination 
of high inlet pressure and low shaft speed, the largest tem- 
perature rises were experienced. 

It can be concluded that any design change that would 
limit or eliminate the source of cooling water must be care- 
fully evaluated as to the thermal implications. Further, in 
any investigation of this nature, due consideration ought to 
be given to the low-speed/hlgh-pressure limiting condition. 
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INTRODUCTION 



Limitation' of the amount of sea water that penetrates 
the hull of a submarine along the propulsion shafting pre- 
sents a variety of problems, most of which have been or 

( 1 2 1 4 ) 

are being examined, ’ s ^’ ’^However, one area of interest, 

in which no experimental work has been done, and little ana- 
lytical effort has been expended, is in the thermal limi- 

to O) 

tations of submarine shaft seals, The type of shaft seal 

presently being utilized in submarine design is of a radial 
face, balanced pressure design with a small amount of leakage 
flow, (see Figure I) A requirement exists to determine the 
quantity of heat generated at the seal interface and the 
efficiency with which this heat can be dissipated. 

Up until quite recently, no heat generation or heat 
dissipation problems were encountered In the design of these 
seals. This was due to the fact that submarines were designed 
with an auxiliary salt water piping system in which the pres- 
sure was maintained slightly above sea pressure. It was only 
a secondary purpose of this system to provide positive flow 
to the shaft seal. This provided a slight amount of leakage 
through the seal and a large amount of bypass flow overboard. 
This latter bypass flow served to effectively remove generated 
heat from the seal. Two conditions have arisen which now make 
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FIGURE I 



SUBMARINE SHAFT SEAL 




Shaft- 



Note: All material 
monel except: 

Seal Ring - Bearium 
Face Ring - Tungsten 
Carbide 
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the thermal limitations of these seals a problem that must be 
faced and quantitatively determined. First, the auxiliary 
salt water system may be completely eliminated in future 
submarine designs. v ' This would effectively remove a 
major portion of the necessary heat sink of the shaft seals. 
Further, as submarine designs and operating profiles trend 
more and more towards deeper submergence, the quantity of 
heat generated at the seal face grows quite rapidly. 

It is the object of this thesis to first present a 
description of an experimental investigation into the ther- 
mal conditions existing in deep submergence submarine shaft 
seals. Then an analytical model of shaft seal is devel- 
oped and utilized to formulate a method of predicting maxi- 
mum interfacial temperatures. In the analysis, some of the 
variables are determined solely from theory, while the de- 
termination of others is based in part on theory and partly 
on the results of the experimental data. The experimental 
results are then presented and compared with the analytical 
results. Extrapolations are made to include conditions be- 
yond the physical capabilities of the test rig. 

It is expected that the analytical procedure determined 
in this paper will be utilized by shaft seal design engineers 
in material and arrangement specifications, and that a great 
deal of costly and time consuming experimental work will 
thereby be avoided. 
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EXPERIMENTAL PROCEDURE 

A 5 5/8" shaft diameter test rig for submarine shaft 
seals was utilized for gathering experimental data. Maxi- 
mum testing conditions were limited to a pressure of 2000 
psi and a rotational speed of 375 rpm. 

A bearium seal ring was selected for instrumentation. 
Bearium (commercial trade name) is a leaded bronze compo- 
sition that has been found quite suitable for use as a 
seal ring in submarine shaft seals, (see Figure II) Five 
1 / 64 ” holes were drilled at equal ( 0 . 044 ") spacings across 
the 0 . 263 " of the seal ring face. These holes were coun- 
ter-sunk and fitted with 1/8" copper tubing for further 
connection to pressure sensing devices, (see Figure III) 
Holes-were drilled for nineteen ( 19 ) chrome alumen thermo- 
couples which were then inserted in the seal ring and seated 
with epoxy, (see Figure IV) 

The seal ring face was then lapped and profiled with 
an X-Y plotter at three separate circumferential positions. 
The face ring (tungsten carbide) and the seal ring were 
inserted in the housing and assembly of the test apparatus 
was completed, (see Figures V and VI) The area ratio (face 
area of the sealing surface to balance area behind face ring 
housing) was 0.523. A sealol step seal, for which leakage 
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FIGURE II 



5-5/8" HIGH PRESSURE DEEP SUBMERGENCE SEAL RING 
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FIGURE V 



SEAL RING INSTALLED IN FRESSURE HOUSING 
( Thermocouple and Pressure Tap Connections ) 
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FIGURE VI 



TEST SEAL AND CONTROL SEAL INSTALLED IN HOUSING 
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( 7 ) 

and torque data was already available from previous testing , w ' 
was inserted in the back half of the test rig. Instrumenta- 
tion was set up to measure the five (5) pressure taps, nine- 
teen (19) thermocouple temperatures, torque, and rpm. (see 
Figure VII) Filtered river water (brackish) was utilized 
as the working fluid. Data was taken with the following 
fixed variables: 



200 


psi 


50 rpm 


0.0 


gpm 


300 


11 


150 " 


0.1 


11 


500 


11 


250 " 


1.0 


11 


1000 


11 


375 " 







1500 " 

2000 « 

Besides the instrumented measurements, inlet pressure, inlet 
water temperature, seal leakage flow, and seal leakage flow 
temperature were measured and recorded. 

Upon completion of all data taking, the seal ring face 
was again profiled at the same locations. The results of 
these profiles (before and after) are shown in Figure VIII. 
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FIGURE VII 



INSTRUMENTATION AND TEST RIG 
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FIGURE VIII 



SEAL FACE PROFILES 
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PRESSURE 

SIDE 



LOW 

PRESSURE 

SIDE 




before test 

after test 
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Ill 



ANALYTICAL PROCEDURE 

Before determining the quantity of heat generated at 
the seal interface and the various parameters which are in- 
volved in that operation, it was first deemed appropriate 
to try to set up an analytical model of the seal apparatus. 
Various techniques were considered and it was finally decided 
that a lumped model based on thermal resistances and a gross 
heat source was best suited for this investigation. Thermal 
resistances were carefully determined based on the most 
accurate predictions of heat coefficients available and 
on judicious estimation of the areas and linear distances 

/ON 

of the relatively complicated seal geometry. ; The con- 
vection quantities were determined by the experimentally 
fixed bypass flow and measured seal leakage flow. (see 
Appendix A) The final model as utilized throughout the 
analysis is displayed in Figure IX. 

The overall heat input at the seal face is a combina- 
tion of a variety of heat source components which vary in 
magnitude with the physical parameters. Basically, there 

are three calculable contributors to the total heat genera- 

( 3 ) 

tion. These are the heat generated by the dissipation 

of pressure energy, the heat generated by viscous shear, 
and the heat associated with boundary lubrication. Thus 
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FIGURE IX 
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the total heat generated at the sealing interface can be 
calculated from the following equation: 



^ ^press + ^visc + ^bound 



or. 



q = 



Q f Ap 



^UAU 

Jc 



fAC Uap 

+ 2 — 



( 2 ) 



In the above equation, only the quantity of flow through 
the sealing faces "Q^" is an experimentally measured quan- 
tity, Pressure drop "^p" and rotational speed "U" are fixed 
variables, while face area "A" is fixed by geometric consider- 
ations and "m" and "J" are constant terms. Therefore, those 
variables remaining to be determined are the coefficient 
estimating the net average contact stress "C " , the seal 
face clearance "c" , and the coefficient of friction " f" . 

The coefficient estimating the net average contact stress 
"C " can be calculated by utilizing the results obtained 

Jr 

from the five (5) pressure taps installed in the seal face. 

The pressure ranges recorded at the five taps spaced at 
radially equal distances apart are plotted in Figures X, 

XI, XII, XIII, XIV, and XV. The range of pressures measured 
for various running conditions of "N" and "Q b " is plotted 
against a non-dimensional constant "K" , where "K" is the ratio 
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FIGURE X 



PRESSURE PROFILE RANGE 
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FIGURE XI 



PRESSURE PROFILE RANGE 
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FIGURE XII 



PRESSURE PROFILE RANGE 
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FIGURE XIII 



PRESSURE PROFILE RANGE 



P 

(psi) 
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FIGURE XIV 



PRESSURE PROFILE RANGE 
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FIGURE XV 



PRESSURE PROFILE RANGE 
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of the difference in radii of the outer seal face and the 
pressure tap to the radial distance between pressure taps. 
Now, since " 0^" is a function of "p", and the pressure 
profiles are known, the seal unbalance may be pictured as: 



P ► 

I | 4 0.263x„ 

1 ;1 0.523 



0.263" 




C 



P 



closing pressure - opening pressure 
inlet pressure 



(3) 



For example, at "p = 2000 psi" , (see Figure X) : 



C_ = 



523p - (p/6) (0.96+0.05+ 0.02+0,01) 

P 



(4) 



= 0.523 - 0.249 = 0 . 2?4 

The balance pressure thus determined and the seal face leakage 
flow are shown plotted against inlet pressure in Figure XVI. 

The seal face clearance "c" for viscous radial flow 
with ”c” and "w" constant can be determined by summing the 
pressure and shear forces across the fluid film and applying 
Newton’s Law of viscous flow. In simple form: 



2 

A£ <tfd u 

w “ dy2 ~ 



(5) 



Since the velocity ”u” is equal to zero at "y=0" and »'y=c". 
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FIGURE XVI 



BALANCE PRESSURE VS. PRESSURE 
and 

SEAL FACE FLOW VS. PRESSURE 
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( 6 ) 




Finally, since 



Q^< = 2ttR 




u dy , 



( 7 ) 




( 8 ) 



In equation (8), "c" was determined by using a weighing 
function on "w" , based on the percent of seal face actually 
experiencing the opening pressure. This weighing function 
was determined from the pressure profiles of Figures X thru 
XV. For "p = 1000 psi", "2/3 w" was assumed; for "p = 200, 
300, and 500 psi", "5/6 w" was assumed; and for "p = 1500 
and 2000 psi", "1/3 w" was assumed. Seal face clearance 
is shown graphically in Figure XVII. 

The determination of a reasonable coefficient of fric- 
tion "f" was finally resolved in the following manner: 



Torque data was available from two sources. Original ex- 
periment produced some data, but difficulty with the me chan 
ics of the rubber coupling at which the torque-meter was 



“ WR " 20.2 C p (psi) 
P 



t_ _ t ( ln-lb) 



(9) 
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6 7 8 9 1 2 3 4 567891 2 3 4 5 67891 2 3 4567891 
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attached made the results somewhat dubious . Data from sim- 

( 7 } 

ilar tests by K. S. Sastelli w provided the author with 
another source of plotted torque evidence. Both of these 
sources were weighed and the combined most accurate predic- 
tion of torque was determined. The coefficients of friction 
thus calculated are shown plotted in Figure XVIII against 
a duty parameter "G" where 

G = 1.224 * 10-9 do, 

P 

With all the variables now accounted for, the total 
heat input at the seal face source was calculated for all 
combinations of the fixed parameters: 

p = 200 - 2000 psi 

N = 150 - 37 5 rpm 

Q b = 0o0 - 1.0 gpm 

From Figures XVI, XVII, and XVIII, the variables were 
extrapolated to include pressures of 3000, 4000, and 5000 psi 
in order to display the effect of very deep submergence on 
the heat input and interfacial temperatures. Total heat 
input, determined as described on the previous pages, is 
shown graphically as a function of inlet pressure and shaft 
speed in Figure XIX. 

With the heat input as source, and Figure IX as model, 
an IBM 7094 computer was utilized to determine the nodal 
temperatures of the resulting 5^5 matrices, (see Appendix B) 
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FIGURE XIX 



HEAT INPUT VS. PRESSURE 
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RESULTS 



The following thirteen (13) graphs make up the body 
of thermal experimental and analytical results of this 
investigation. In Figure XX is shown the actual seal face 
temperatures (determined analytically) versus the source 
heat input. Figures XXI and XXII show actual seal face 
temperature (determined analytically) plotted against inlet 
pressure for bypass flows of 0.0 and 0.1 gpm respectively. 
Figures XXIII through XXXII depict both an analytical and 
an experimental non-dimensional temperature quantity "T" , 
where 



T = 




( 11 ) 



This is plotted against "Q^" , the bypass flow rate, for 
varying inlet pressure from 200 to 5000 psi. Note that 
Figure XXVIII is a duplicate presentation of Figure XXIX 
with the exception of the scale utilized. Note also that 
all data in Figures XXX, XXXI, and XXXII is from the analyt- 
ical extrapolation. 
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FIGURE XX 



SEAL FACE TEMPERATURE VS. HEAT INPUT 
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FIGURE XXI 



SEAL FACE TEMPERATURE VS. PRESSURE 
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FIGURE XXII 



SEAL FACE TEMPERATURE VS. PRESSURE 
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DISCUSSION OF RESULTS 



The first three figures of Chapter IV present, in cap- 
sule form, the major points learned, from this investigation. 
The most obvious and also most important result is best seen 
in Figure XX on page 30. With no bypass cooling water flow, 
the seal interfacial temperature rises sharply with heat in- 
put, so swiftly in fact that a heat input of only about 
1600 BTU/HR is required for vaporization of the film. How- 
ever, it can be seen that only a very slight quantity of 
bypass cooling water flow (0.1 gpm) will decrease the tem- 
perature rise to a safe level for even the highest predicted 
value of heat input. Further, if the quantity of bypass 
cooling water flow is increased to 1.0 gpm, the seal inter- 
facial temperature rise is almost negligible for even large 
increases in heat input. 

Figures XXI and XXII on pages 33 and 34 respectively 
display the same phenomena when examined together. Taken 
separately, some further discussion presents itself. In 
Figure XXI, it can be seen that for zero bypass cooling water 
flow, the limiting seal inlet pressure is between 3400 psi 
and 3700 psi, depending on the rotational speed of the ap- 
paratus. A very interesting phenomenon also discernable in 
this Figure is the sharper rise in temperature of the low 
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speed analytical data at high values of input pressure. 

This can be attributed first to the results of Figure XIX 
on page 28, which in turn is really a function of the three 
heat input variable terms. Further examination of the ana- 
lytical data shows that the key lies in the heat associated 
with boundary lubrication, the third term in the heat equa- 
tion (Equation (2) on page 15). For example, the low value 

-10 

of duty parameter "G=l. 34x10 " at "p=5000psl" and "N=150rpm" 

results in a relatively high value of the coefficient of 
friction "f=0„0l8", as can be seen on the extrapolated por- 
tion of Figure XVIII on page 27. 

The same phenomenon is also discernable (to a lesser 
extent) in Figure XXII on page 32. Here, with a small amount 
of bypass cooling water flow, the temperature rises are not 
very sharp, and are also not significantly dependent on RFM. 

At very low pressures, the spread of results between 
the experimental and the analytical data was a maximum, (see 
Figures XXIII, XXIV, and XXV) However, at higher values of 
the seal inlet pressure (p = 1000 psi and above) the experi- 
mental and analytical results show a marked similarity in 
trends and values. It is extremely interesting to note that 
the trend for the combination of low RPM and high pressure 
to result in a greater heat input and hence higher inter- 
facial temperature begins at a considerably lower pressure 
in the experimental results then it does in the analytical data. 
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The fact that this trend did exist both experimentally and 
analytically is quite significant. 

The curves of Figure XXVIII on page 38 display in true 
focus the sharp rise in the temperature associated with de- 
creased bypass cooling water flow. It would have been de- 
sirable to display all of the temperature increases in this 
form, but this would have resulted in the curves being so 
bunched as to be unreadable; hence the logarithmic displays. 



45 



VI 



CONCLUSIONS AND RECOMMENDATIONS 

Deep submergence submarine shaft seals of a similar 
design to those presently being used require a minimum 
amount of bypass cooling water flow. Any design change 
that would limit or eliminate the source of cooling water 
must be carefully evaluated as to the thermal implications 
involved. It is recommended that the method of evaluation 
should follow that presented in the Analytical Procedure 
and in the Appendices of this paper. 

Further, in any investigation of this type, due consid- 
eration should be given the low-speed/high-pressure limiting 
condition. This warning is especially emphasized since one 
would intuitively expect that a high-speed/high-pressure 
limiting condition would be all that required examination. 

In future experimental investigations with submarine 
shaft seals, it is recommended that a more satisfactory method 
of gathering accurate torque data be devised and utilized. 

Only then can a great deal of confidence be given the computed 
coefficients of friction and hence heat source inputs. 

Finally, computer programs similar to that described in 
Appendix B are available in most large computation centers, 
and ought to be utilized in investigations of this nature. 
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APPENDIX A 



DETERMINATION OF THERM AI RESISTANCES 

The thermal resistances shown diagrammatlcally in Fig 
ure IX on page 14 were calculated from the following basic 
equations : 



q = h A &T 

k A aT 

q « — 

q = Q C aT 



R = 



h A 



R = 



R = 



k A 

1 

Q C 



(12) 

(13) 
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The following constant coefficients were established: 





1: BTU 


h btu 


c PT y 




- hr ft °F 


- hFTt^oF 


- lb 6 F 


Bearium 


20.4 






Monel 


15.0 






Tung . Carb . 


30.0 






Epoxy- 


0.6933 






Cooling Water 


0.33 




1.0 


Monel-air 




0.8 




Interface 








Bearium-Monel 




1000.0 




Interface 








Monel-Water 




3570.0 




Interface 








Bearium-Water 




5950.0 




Interface 









The value of "Q^" was set either at 0*0 (closed cir- 
cuit) or at Ool or 1.0 gpm. Values of "Q f " were determined 
individually for the range of pressures utilized, (see Fig- 
ure XVI on page 23) 

The following cross-sectional areas of heat transfer 
surface and lengths of heat transfer paths were estimated 
by examination of the particular geometry of the seal under 
study. A similar analysis would be required for any modi- 
fication In seal geometry. However, for a simple scaling 
modification, the areas and linear distances could be simi- 
larly scaled. 





A (in 2 ) j 


L (in) | 


Eearium 


31.0 


1 . 5 | 


Bearium-Monel Interface 


71.0 




Monel (seal ring housing) 


71.0 


1.2 


Monel-Air Interface 


318.0 




(seal ring housing) 
Bearium-Water Interface 


34.0 




Sealing Water at Face 


2.1 


1.4 


Sealing Water at Leakoff 


1.0 


1.8 


Seal Water-Air Interface 


CO 




Tungsten Carbide 


13c0 


0.5 


Epoxy 


34.0 


0.005 


Monel (face ring housing) 


| 22.0 
! 


1.0 


Monel (shaft) 


| 150.0 


O 

o 

o 
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A (in 2 ) 


L (in) 


Monel-Air Interface 


oo 




( shaft) 

Monel-Water Interface 


60 . 0 




(face ring housing) 
Bypass Water at Face 


7.0 


1.6 


Ring Housing 
Water-Monel Interface 


157.0 




(seal housing) 
Monel (seal housing) 


171.0 


2.0 


Monel-Air Interface 


472,0 




(seal housing) 
Bypass Water at 


14.0 


10.0 


Bypass Channel 
Bypass Water-Air Interface 


<=*=> 
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APPENDIX B 



COMPUTATIONAL PROCEDURE 



The computer program utilized in this thesis can be 
found in the Department of Naval Architecture and Marine 



Engineering Computer Library at the Massachusetts Institute 
of Technology under the coding of NA 021,0. The computer 



listing 


follows : 




FORTRAN 




DIMENSION A( 1 1 , 1 1 ) ,B( 1 1 , 1 1 ) ,X( 1 1 ) , E( 11 ) , IL( 20 ) 


500 


READ 400, JOB 


400 


FORMAT (15) 




IF (JOB) 501, 502, 503 


501 


CALL EXITM 


502 


CALL EXIT 


503 


IF ( JOB-2) 100,200, 300 


c 


FART 1 SIMULTANEOUS EQUATION SOLN BY XSIMEQF 



C F FACTOR BY WHICH XSIMEQF MULT. VALUE OF DET. A(I,J) 
C JJ-NUMBER OF COLUMNS IN MATRIX B 

C M=1 SOLN GOOD, BRIGHT BOY.A(I.J) REPLACED BY X(I,J) 

C M=2 OVERFLOW, BROKE THE BANK 
C M=3 MATRIX A IS SINGULAR, GO TO JAIL 
100 READ 15 
PRINT 15 

READ 10, N, ( (A(I,J) , 1=1, N) , J=1,N) 

PRINT 1 1 ,N 
DO 16 1=1, N 

PRINT 19, ( A( I , J) , J=1,N) 

16 CONTINUE 

PRINT 24 

READ 2, (B(K, 1) ,K=1,N) 

PRINT 2, (B(K, 1) , K= 1 , N ) 

2 FORMAT (5F13.4) 

PRINT 1? 

PRINT 18 
PRINT 24 

10 FORMAT (I5/(5F13.4) ) 

11 FORMAT ( 4H N= ,13) 

12 FORMAT ( 16H GOOD BOY DET = , 1PE14. 6/( op5Fl4. 5) ) 

13 FORMAT (16H BROKE THE BANK) 
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14 FORK AT (54H GO TO JAIL, DO NOT PASS GO, DO NOT COL- 

LECT 200 DOLLARS) 

15 FORMAT (54H ) 

17 FORMAT ( 1H2) 

18 FORMAT ( 8H ANSWERS) 

19 FORMAT ( 1H0.10F11.4) 

24 FORMAT (1H0) 

25 FORMAT (F10.4) 

26 FORMAT (8H SCALE= ,F10.4) 

READ 25, F 

PRINT 26, F 
PRINT 24 
JJ=1 

M=XSIMEQF( 11,N,JJ,A,B,F,E, ) 

GO TO (20,21,22) ,M 

20 PRINT 12, F, (A(I, 1) ,1=1, N) 

GO TO 23 

21 PRINT 13 

GO TO 23 

22 PRINT 14 

23 GO TO 500 



The program, when supplied with the proper input control 
cards, will solve simultaneous algebraic equations up to a 
maximum of 10 equations and 10 unknowns: 

a. When the solution is successful - BRIGHT BOY, 

b. When overflow occurs - BROKE THE BANK, 

c. When the matrix is singular - GO TO JAIL. 
d 0 When matrix to be inverted is not square - 

(caustic comment). 

A complete write-up of this FAP program can be obtained from 
the MIT Consultants (cc-174-3). 

INPUT 

Card 1 Format (15) 

This is the control card. Place a 1 In column 5. 

Card 2 Format ( 54H )_ 

Flace a 1 In column 1 of the data card, then type 
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any comment that you wish to have printed. 



Card 3 Format (15) 

The number of unknowns in the set of simultaneous 
equations is punched in columns 4 and 5i the units 
digit appearing in column 5. 

Group 4 Format (5F13.4) Cards = ' ■ ^ ' 

The value of the coefficients are punched five to 
a card. The coefficients are punched column by col- 
umn (i.e. , first column - all rows, second colume - 
all rows, etc.). 

Group 5 Format (5F13.4) Cards = ^ 

The value of the algebraic equations are punched 
five per card. There is no need to fill all 5 
fields on the card. 

Card 6 Format (FI 0.4) 

This value serves as a scale factor by which XSIMEQF 
multiplies the value of the determinent of the mat- 
rix. If the value is one, a 1.0 is punched anywhere 
in the above field. The scale factor is printed out 
in the output data as SCALE = . This scaling option 
applies to those who are manipulating large numbers. 



Sample input and output data follow. The comment card 
(1st card) indicates in order: inlet pressure (psl), bypass 

flow (gpm), and rotational speed (rpm). 
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200 1.0 375 
N= 5 



576.40 


-35.00 


-10.00 


-18.40 


-513.00 


-35.00 


1436.72 


-1400.00 


0 . 


0 . 


- 10.00 


-1400.00 


1420.00 


0 . 


0 . 


-18.40 


0 . 


0 . 


1537.10 


-1500.00 


-513.00 


0 . 


0 . 


-1500.00 


2529.00 


505.20 


148.00 


860.00 


1609.00 


44400.00 



ANSWERS 
SCALE= IcOO 

GOOD BOY DET = 3.619504E 13 

87.86 87.48 87.47 86.95 86.95 



1000 0.1 150 

N= 5 



110.85 


-35.00 


-6.05 


-18. 40 


-51.40 


-35.00 


1436.72 


-1400.00 


0 . 


0 . 


-6.05 


-1400.00 


1412.05 


0 . 


0 . 


-18.40 


0 . 


0 . 


1537.10 


-1500.00 


-51.40 


0. 


0 . 


-1500.00 


1605.20 


285.70 


148.00 


516.00 


1609.00 


4629.99 



ANSWERS 
SCALE= 1.00 

GOOD BOY DET = 6.280575E 11 

92.84 91.77 91.75 89.38 89.38 



1 
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3000 0 250 
N= 5 



54.24 


-35.00 


-0.72 


-18.40 


-0.12 


-35.00 


1436.72 


-1400.00 


0 . 


0 . 


-0.72 


-1400.00 


1401.41 


0 . 


0 . 


-18.40 


0 . 


0 . 


1537.10 


-1500.00 


-0.12 


0 . 


0 . 


-1500.00 


1502.72 


1067.30 


148.00 


59.30 


1609.00 


224.00 



ANSWERS 
SCALE= 1 . 00 

GOOD BOY DET = 3.888511E 10 

173»8o 168.25 168 0 21 126.87 126.80 



In all, nlnty-one (91) data points (solution of 91 
separate 5^5 matrices) were obtained utilizing this pro- 
cedure. The five temperatures given in the output of 
this program conform respectively to the five lettered 
(A,B,C,D,E) nodes in Figure IX on page 14. 
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